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Abstract 

With the BABAR detector at the PEP-II asymmetric B Factory at SLAC, we present the first 
search for the decays B° — » £ + l~j (£ = e, fi). Using a data set of 292 fb -1 collected at the T(4S) 
resonance, we find no significant signal and set the following branching fraction upper limits at 
90% confidence level: B(B° -> e + e~j) < 0.7 x 10~ 7 and B(B° -> /i + /i"7) < 3.4 x lO" 7 . 
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1 INTRODUCTION 



Rare decays induced by flavor changing neutral currents occur at loop level in the Standard Model 
(SM) and they are sensitive to the flavor structure of the SM as well as to the new physics beyond the 
SM. Studying radiative B meson decays such as B° — > £ + £~j 5 can provide us essential information 
on the parameters of the SM, such as the elements of the Cabibbo-Kobayashi-Maskawa matrix and 
the leptonic decay constants. As explained below, the branching ratios that the SM predicts are 
much below the experimental sensitivity. Hence this paper presents a search for new physics. 

The most important contribution to B° — > £ + £~j arises from radiative corrections to the pure 
leptonic processes B° — > £ + £~ which suffers from the helicity suppression. The short distance 
contributions to B° — > £ + £~ come from the box, Z-boson and photon-mediated diagrams. If a 
photon is emitted from the final charged lepton lines, the amplitude is proportional to the lepton 
mass me and is helicity suppressed. When a photon is attached to any charged internal line, 
the contributions are also significantly suppressed by a factor of m^/rn^. Therefore, the main 
contribution is when a photon is radiated from the initial quark lines as shown in Fig. ^ 




(a) (b) (c) 

Figure 1: Feynman diagrams for the main contributions to B° — > £ + £~*y decay. The signal photon 
can also be radiated from the initial d quark line. 

The expected branching fractions are summarized in TableHfor both radiative and non-radiative 
dilepton B decays P^. This note presents the first search for B° — ► £ + £~^/. 

Table 1: The predicted Standard Model branching fractions. 



Channel 


Branching Fraction 


B l) -» e + e"7 
B° -> 


1.5 - 4 x l(T iU 
1.2 - 3 x 1(T 10 


B u -> e + e- 
B° -» fi + fi- 


2.1 x l(T i!3 
9 x l(T n 



Throughout the document, I denotes either e or fi, and the charge conjugated states are implicitly included. 
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2 THE BABAR DETECTOR AND DATASET 



The data used in this analysis were collected with the BABAR detector at the PEP-II storage ring 
at the Stanford Linear Accelerator Center and correspond to an integrated luminosity of 292 fb _1 
accumulated at the T(4S) resonance, which is equivalent to 320 million BB events, and 27 fb -1 
accumulated at a center-of-mass (CM) energy about 40MeV below the Y(4«S) resonance. 

The BABAR detector is described elsewhere The 1.5 T superconducting solenoidal magnet 
contains a charged particle tracking system with a silicon vertex tracker (SVT) followed by a 
drift chamber (DCH), a ring imaging Cherenkov detector (DIRC) dedicated to charged-particle 
identification, and an electromagnetic CsI(Tl) calorimeter (EMC). The segmented flux return (IFR) 
is instrumented with resistive plate chambers. About 1/3 of these chambers have been replaced 
with limited streamer tubes which provide higher muon identification efficiency. This change affects 
the most recent 77 fb _1 of data. 

A full BABAR detector Monte Carlo (MC) simulation based on GEANT4[Hj is used to evaluate 
signal efficiencies and to identify and study background sources. 



We reconstruct B° candidates with two oppositely-charged leptons and a photon. The leptons are 
required to originate from a common vertex, and the B° candidate is required to be consistent with 
coming from the beam interaction point. Since the signal events contain two neutral B mesons and 
no additional particles, the total energy of each B meson in the CM frame must be equal to half 
of the total beam energy in the CM frame. We define 



where E£ is the (e + or e~) beam energy in the CM frame, p? and mi are the momentum in the 
CM frame and the mass of the daughter particle i (i = £ + ,£~ ,7), respectively. In Eq. ^ -^beam * s 
used instead of the B meson energy in the CM frame because E^ eam is more precisely known. For 
correctly reconstructed B° mesons, tues has a maximum at the nominal B° mass with a resolution 
of about 3 MeV/c 2 and AE is near zero with a resolution of about 30MeV. 

The B° -» £ + £--y candidates are selected in the -0.5 < AE < 0.5 GeV and 5.2 < m E s < 
5.3GeV/c 2 range. The size of the signal box is chosen to be approximately ±3cr of the AE and 
m E s distributions: \AE\ < 0.123 GeV and 5.27 < m E s < 5.288 GeV/ c 2 . The resolutions in AE 
and nriES are obtained from fits to the signal MC distributions assuming Crystal Ball function^ 
shapes. A larger region which contains the signal box was blinded during the development of this 
analysis: \AE\ < 0.164 GeV and 5.267 < rriEs < 5.29 GeV/ c 2 . The remaining region, the sideband 
area, is retained for background studies and data/MC comparison. 

To minimize the number of mis-identified particles, the leptons are required to satisfy stringent 
electron and muon identification criteria. Electrons are identified using a likelihood method with 
EMC and DIRC information. The electron identification efficiency is about 93% and a pion fake 
rate is less than 0.1%. Muons are identified using a neural network method with IFR information. 



3 ANALYSIS METHOD 




(i) 



(2) 
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The muon identification efficiency is about 70% and a pion fake rate is less than 3%. For photons, 
we require the shower lateral moment 5 to be less than 0.53(0.51) for the electron(muon) mode. 

Leptons and photons are required to respect strict acceptance criteria. We have photon quality 
requirements on the number of crystals(> 10) and photon energy(> 0.3 GeV) for the electron mode 
only. These requirements help reject beam background, the initial state radiations, and higher 
order QED backgrounds which are not modeled in MC for the electron mode. 

We apply vetoes on any leptons which may come from J/ifi or ip(2S) decays. These vetoes are 
chosen in the 2-dimensional plane of AE and the invariant mass of the two leptons, m^L- The J/ip 
veto region is the union of the following three regions in the AE-m^i plane for electron(muon) 
mode: 

• 2.90(3.00) < m DL < 3.20 GeV/ c 2 , 

• for moL > 3.20 GeV/ c 2 region, a band in the AE — tudl plane defined by 
1.11c 2 x m DL - 3.58(3.53) GeV < AE < 1.11c 2 x m DL - 3.25(3.31) GeV, 

• for rri£)L < 2.90(3.00) GeV/c 2 region, a triangle in the AE — tudl plane defined by 
AE < 1.11c 2 x m DL - 3.25(3.31) GeV. 

Photon candidates whose invariant mass with other photons in the event are in the range of 0.115 < 
m 77 < 0.155 GeV/ c 2 are vetoed to remove photons from ir° decays. We studied the effect of an rj 
veto and found it to be negligible. 

We require that the invariant mass of dilepton system is between 0.3 and 4.9(4.7) GeV/c 2 for 
the electron(muon) mode to reject non-BB background. Further suppression of background from 
non-BB background is provided by a series of topological requirements. We require | cos 9t\ < 0.8, 
where 6t is the angle in the CM frame between the thrust axis of the particles that form the 
reconstructed B° candidate and the thrust axis of the remaining tracks and neutral clusters in 
the event. We have a requirement on the ratio of the second to zeroth Fox- Wolfram moment [0] of 
i?2 < 0.31. We define a Fisher discriminant based on the following variables: the angle between 
the B° direction and the beam axis, the angle between the thrust axis of the B° candidate and the 
beam axis, and the summed momentum of the rest of the event tracks and neutrals in nine non 
overlapping volumes delimited by cones centered around the thrust axis of the B° candidate with 
half-angles in 10° steps. The distributions of the Fisher discriminants are shown in Fig. |5J The 
arrows show the allowed regions. If multiple B° candidates pass all the selection requirements in 
an event, only the one with AE closest to zero is retained. The average number of B° candidates 
per event is 1 (1.05) for the electron (muon) mode. 

To assess potential background contributions peaking like the signal in AE and m^s, we exam- 
ined 32 exclusive hadronic and semileptonic decay modes, and found no significant contribution. 

The requirement levels are optimized by minimizing the expected upper limit branching frac- 
tion at 90% confidence level (C.L.). The upper limit is calculated using the Feldman-Cousins 
method 0. The normalizations and shapes of the data and MC distributions are compared at 
various requirement levels in the sideband area. They show reasonable agreement at all levels. 

To estimate the background level in the signal box, three different sideband boxes are used, as 
indicated in Table El We perform an unbinned maximum likelihood fit on the tues distribution in 
the combined upper and lower sideband boxes with an Argus function jSj., as shown in Fig. El The 
end-point of the Argus function is fixed at 5.29 GeV/ c 2 . We use this parameterization to extrapolate 
the background level found in the middle sideband (12 events for B° — > e + e~7 sample and 13 events 
for B — > fi + fj,~j sample) into the signal box. The error on the expected number of background in 
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Figure 2: The distribution of the Fisher discriminant for selected events in the signal box for (a) 
B° — > e + e~7 and (b) B° — > / u + //~7 sample. The points show the predicted distributions for signal 
events, the histograms show the various background contributions. The J/ifi , ip(2S) and 7r° vetoes 
are applied. The signal distribution is normalized to the branching fraction of 1 x 10~ 6 . 



the signal box from the fit is estimated by varying the Argus parameters by ±lcr. The background 
has been estimated with alternative methods, using sidebands in AE, and repeating the fit after 
relaxing the selection requirements. The estimates obtained are in general compatible within errors. 
The fit to the tues sidebands was chosen as the preferred method because it yields the smallest 
statistical uncertainties. 

Table 2: Definitions of the different sideband boxes used. 



Sideband Box 


span in AE ( GeV) 


span in txies (GeV/c 2 ) 


Upper Sideband 
Lower Sideband 
Middle Sideband 


(0.164, 0.5) 
(-0.5, -0.164) 
(-0.123, 0.123) 


(5.2, 5.3) 
(5.2, 5.3) 
(5.2, 5.26) 



4 SYSTEMATIC STUDIES 

Systematic uncertainties are evaluated using independent data control samples. The dominant sys- 
tematic uncertainty is related to the reconstruction of the signal photon energy which is determined 
using e + e~ — > /i + /i~7 decays. The uncertainty is 1.6% for both modes. The systematic uncertainty 
from the particle identification has been determined using an independent control sample of J/i/j 
decays. It is 0.7% and 1.3% for the electron and muon mode respectively. The uncertainty on the 
number of BB events is The tracking efficiency is determined from e + e~ — > t + t~ interac- 
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Figure 3: An unbinned maximum likelihood fit on uies distribution in the combined upper and 
lower sideband boxes, using an Argus function for (a) B° — > e + e~7 and (b) B° — > ^x + /i~7 sample. 
The points represent data and the solid lines show the fit function. 



tions, with one tau decaying leptonically and the other to three charged hadrons. The uncertainty 
is 0.94% for both electron and muon modes. All contributions to the systematics are added in 
quadrature to give a total relative systematic uncertainty of 2.3% for the electron mode and 2.5% 
for the muon mode on the branching fraction. 



5 RESULTS 

As shown in Fig. 0] and Table El zero and three events were found in the electron and muon 
modes, respectively. The number of events found in the signal box is compatible with the expected 
background for both modes. 

An upper limit on the branching fraction is computed using 



N B o ■ e. 



(3) 



siy 



where Nul is the 90% C.L. upper limit for the signal yield as determined by the method described 
in ^U] with both statistical and systematic errors, N B o is the number of the neutral B mesons and 
€ s ig is the signal reconstruction efficiency from the signal MC sample. N B o is equal to the number of 
BB events, as we are assuming B(T(45) -> B°B°) = B(T(45) -» B+B~). The obtained 90% C.L. 
branching fraction upper limits are B(B° — > e + e~7) < 0.7 x 10 _7 and B(B° — > fi + < 3.4 x 10 -7 . 



6 CONCLUSION 



A search for the B° — 
BB events. We obtain 



= e or fi) modes has been performed based on 320 million 
o C.L. upper limits for the branching fraction of B(B° — ► e + e~j) < 
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Table 3: Summary of the results where n f, s and n^ p are the observed and expected number of 
background events in the signal box, e s i g is the efficiency, Njjl is the 90% C.L. upper limit for the 
signal yield, and Bul(B — > t+l"^) is the upper limit on the branching fraction at the 90% C.L. 
Systematic uncertainties on n e ^ p and e s i g are given. 



Decay Mode 


Tlobs 


exp 


e« 9 (%) 


N UL 


b ul (b» -> e+e-y) 


e + e~ 



3 


1.28 ±0.80 
1.40 ± 0.42 


6.07 ±0.14 
4.93 ±0.12 


1.33 
5.30 


0.7 x 10- y 
3.4 x 10~ 7 



> 

03 

O0.4 
LU 

< 

0.2 




• 


BaBar 

preliminary 


> 

03 

O0.4 
LU 

< 

0.2 


• 

— • 

• 


BaBar 

preliminary 

« 
















• ' 











* • • 
• • 









* . • 

_• • * 
• 


• 








■ — — 


• * 












T« * 

♦ • 






-0.2 


• • 


• • 








-0.2 




• 


9 




-0.4 


• • 

- 


• 

• 

<• 


» • 






-0.4 


.. • ' 


/ 

• • A* . 

A 







5.2 5.22 5.24 5.26 5.28 5.3 5.2 5.22 5.24 5.26 5.28 5.3 

m ES (GeV/c 2 ) m ES (GeV/c 2 ) 

(a) (b) 

Figure 4: Distribution of events in tues and AE for (a) B° — > e + e~7 and (b) i?° — > / Lt + / Lt _ 7 sample. 
The solid box represents the signal box. The three dashed-lined area show the upper, middle, and 
lower sideband boxes from top to bottom. 



0.7 x 10 7 and B(B° — > 7) < 3.4 x 10 7 . These are the only limits currently available for 

these decay modes. 
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